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Abstract 

The radiative decay of the next-to-hghtest neutrahno into a lightest neutrahno and a photon is 
analyzed in the minimal supersymmetric standard model. We find that significant regions of the 
supersymmetric parameter space with large radiative branching ratios (up to about 100%) do ex- 
ist. The radiative channel turns out to be enhanced when the neutralino tree- level decays are 
suppressed either "kinematically" or "dynamically". In general, in the regions allowed by data 
from CERN LEP and not characterized by asymptotic values of the supersymmetric parameters, 
the radiative enhancement requires tan/3 ~ 1 and/or Mi ~ M2, and negative values of fi. We 
present typical specific scenarios where these necessary conditions are satisfied, relaxing the usual 
relation Mi = | tan'^0wM2. The influence of varying the top-squark masses and mixing angle when 
the radiative decay is enhanced is also considered. Some phenomenological consequences of the 
above picture are discussed. 
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1 Introduction 



Neutralinos {Xi, i = 1,...4; m^o < ... < tri^o) are among the lightest supersymmetric 
(SuSy) partners predicted in the Minimal Supersymmetric Standard Model (MSSM) [|I|].[] 
In particular, the Xi is usually the lightest SuSy particle (LSP). Hence, for conserved im- 
parity, Xi is always present among the decay products of any superpartner, giving rise to 
large amounts of missing energy and momentum in the final states corresponding to pair 
production of SuSy particles in e^e~ /pp/ep collisions. The heavier neutralinos have in 
general a rather complicated decay pattern towards the LSP, with possible intermediate 
steps involving other neutralinos and/or charginos as well as two-body decays with on-shell 
Z° /W^^ Higgses or sfermions in the final states Being the lightest visible neutralino, 
the next-to-lightest neutralino X2 is of particular practical interest. It would be among the 
first SuSy partners to be produced at the CERN e^e~ collider LEP2 and at the Fermilab 
TeVatron 0, |. 

The dominant X2 decay channels are, in general, tree-level decays into a lightest neutralino 
and two standard fermions through either a (possibly on-shell) Z° or sfermion exchange. 
Accordingly, the X2 decays into a Xi plus a vivi or qq pair. When X2 is heavier than 

the lightest chargino, also cascade decays through a Xi can become relevant. Hence, one can 
have, as final states, t^i'~vi/pgi^ i'^u^qq', qiq[q2<l2, plus a Xi- An additional possibility is the 
two-body mode xl ~^ Xih°{A°), when the Higgs boson(s) is (are) light enough |^. This 
may give rise to events with a bb/r^r^ pair and missing energy and momentum (^ and ]/) 
in the detectors. 

The radiative X2 decay into a photon and a lightest neutralino X2 ~^ Xil provides a 
further decay channel with an even more interesting signature. Analytical formulas for the 
corresponding width can be found in Ref.0.0 Because of the higher-order coupling, this 
channel is characterized in general by rather low branching ratios (BR's). Nevertheless, a 
comprehensive study of the X2 decay rates shows that there are regions of the SUSY 
parameter space where the radiative decay is important and can even become the dominant 
X2 decay. 

In this paper, we perform a detailed phenomenological analysis of the decay X2 ~^ X1I 
in the MSSM. In particular, we analyze the regions of the SuSy parameter space where 
the radiative decay is enhanced. Some asymptotic regimes that give rise to large radiative 
BR's have been considered in Ref.0. In these particular considerable hierarchy 

is present among the different mass parameters of the MSSM Lagrangian, that by now 
corresponds to regions of the SuSy parameter space either partly or entirely excluded by 
LEP searches. In our study, we go beyond such asymptotic scenarios and consider regions of 
the SuSy parameter space where the parameters Mi 2 and \fi\ have values roughly included 
in the range [Mz/4:,4:Mz]- We relax the usual condition on the electroweak gaugino masses 

^In this paper, by MSSM we refer to the supersymmetric extension of the Standard Model with minimal 
gauge group and particle content. No additional assumptions (such as unification assumptions at a large 
scale) are understood, unless exphcitly stated. 

^In checking the results of Ref. we found that the first line in Eq.(59), p. 281, should be more properly 
written as follows: 

Gl = 2coset[T3tZl-^ + ettane,,.Za] + sm9tZ,4 

mw sm p 

In this way, it can also be correctly generalized to the f ~ b,T, . . . case, by substituting T^t, et,Qt ~^ '^sfT^fj^f- 
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Ml = |tan^^vi/M2, which holds [through one-loop renormalization group equations (RGE's)] 
when one assumes their unification at a scale Mqut ~ 10^^ GeV, where the gauge couplings 
assume the same valuej^ In this paper, we will treat Mi 2, as well as the masses of the 
individual sfermions, as low-energy, free parameters. As for the sfermion masses, we will 
assume, for simplicity, some common value (or two different values for sleptons and squarks), 
whenever the individual mass values are not particularly relevant for the analysis. 

As we will show, taking both Mi and M2 as free parameters can produce new interesting 
scenarios beyond the ones already considered in Ref. 0. In that paper, we already singled 
out some regions of the parameter space where the tree-level X2 decays are suppressed and 
BR(x2 x\l) is large. There, we partly misinterpreted the origin of the suppression in a 
few points of the most promising regions in the (/i, M2) plane, ascribing it to the particular 
physical composition of Xi2- ^ill be clear from the following more complete analysis 
that in those scenarios a kinematical suppression of the tree-level decays can be effective, in 
addition to a dynamical one. We indeed noticed in Ref. 0, that the two hghtest neutrahnos 
are nearly degenerate in the interesting cases and that their mass difference {m^o—rri-^o) grows 
monotonically with tan/5. Now, we will show how different mechanisms can contribute to 
the radiative BR enhancement. These scenarios will be thoroughly analyzed by a systematic 
investigation of the regions of the SuSy parameter space where a large BR(x2 ^ Xi7) 
regime may be present. For instance, a dynamical suppression of the x\ tree-level decays 
occurs when Xi and X2 have a different dominant physical composition in terms of gauginos 
and Higgsinos. 

The latter scenario can have particular relevance for explaining events like the e+e"77-|- 
event recently observed by the Collider Detector at Fermilab (CDF) Collaboration at the 
Tevatron 0. That event, characterized by the presence of hard photons, electrons and large 
missing transverse energy, can be hardly explained within the Standard Model (SM). A pos- 
sible solution to this puzzle can be found within the minimal SuSy models, by interpreting 
the CDF event as a result of selectron- (or chargino-) pair production, provided a large BR 
for the neutralino decay into a LSP and a photon is predicted. More generally, the presence 
of a large radiative neutralino decay BR is crucial to obtain, within SuSy high rates for 
final states associated to a signature of the kind 77 -|- X+ ^, where X = t^l^'^~ , qq, ... or 
nothing. 

Two different SuSy "models" have been proposed up to now to explain the e+e"77+ I^t 
CDF event [l^l-S The first one arises within theories with low-energy supersymmetry 
breaking, where the breaking is transmitted to the visible sector by nongravitational interac- 
tions [0. In such a scenario, the gravitino G turns out to be naturally the LSP and, if light 
enough (i.e., for ^ 1 keV), it can couple to standard SuSy particles strongly enough to 
be of relevance for collider phenomenology |13[. If one assumes that the lightest standard 
SuSy particle is still the lightest neutralino and -R-parity is conserved, all the heavier SuSy 
particles produce, at the end of their decay chain, (at least) one Xi- Then, Xi decays ra- 



■^Note that this not necessarily corresponds to relaxing all the gaugino mass unification conditions. One 
can still imagine that M2 and M3, that is the parameters which correspond to non-abelian gauge groups, 
unify in the usual way, while the unification relation between Mi and AI2 may be different from the usual 
one and unknown. 

^As we were completing this paper, three other papers appeared which discuss the CDF event in various 
contexts pl| . 
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diatively, Xi ^ G^y, through its photino component, if the latter is not tuned to zero 
Assuming the Xi radiative decay occurs well into the detector (i.e., close enough to the main 
vertex, which requires m^, ^ 250 eV), such a model can provide a satisfactory explanation of 
the e'''e"77+ event ||, ||. Once low-scale SuSy breaking and very light gravitino scenar- 
ios are trusted, the presence of a large-BR xl radiative decay and the consequent signature 



of hard and central photons and missing energy are almost automatic [§, ^ |T^. Although 



quite general, such a hypothesis does not allow to predict much about the SuSy parameters 
apart from rriQ. In particular, an interpretation of the CDF event within this framework can 
single out some ranges for the physical masses of the involved particles only on the basis of a 
careful analysis of the kinematical characteristics of the event |T0[. No specific statements 
about the values of the parameters in the SuSy Lagrangian (Mi, M2, /i, tan/5, etc.) and, 
hence, no detailed predictions of the general related collider phenomenology can be achieved. 

Somehow opposite is the situation if the CDF event (or a general 77 -|- X-\- ^ event) 
is explained within the MSSM, where the gravitino is heavy, the lightest neutralino is the 
LSP and the hard photons and the missing energy are due to the one-loop xl ~^ Xi7 decay 
H, |TOl|. In this certain adjustment of the MSSM parameters is required (both in 

the dynamical- and in the kinematical-enhancement scenarios) in order to get large radiative 
BR's and large rates for events with hard central photons and missing energy. Also, if the 
hard photons are emitted by rather soft xV^^ ^^e (i?/namzca/ enhancement is the most effective 
mechanism in this respect. Then, one is in general able to select rather narrow ranges of the 
SuSy parameters, if the 77 + X+ ^ events are interpreted in this framework [||, [l^. Hence, 
such a framework can be quite predictive also about the SuSy collider phenomenology that 
should show up in the future. 

In this work, we intend to investigate the latter hypothesis in a general framework. We 
simply look for regions in the usual MSSM parameter space where the X2 ~^ Xi7 decay has 
sizeable BR's. We find that in order to have a large BR(x2 Xil) (up to about 100%) 
one needs in general tan/5 ~ 1 and/or Mi ~ M2, in addition to fi < 0.Q This is a quite 
general requirement, while further conditions on the gaugino mass parameters Mi, M2, the 
Higgsino mass \fi\, and tan/? can guarantee either a dynamical or a A;mematzca/ enhancement 
of BR(x2 Xil)- These two possibilities will give rise to rather different spectra for the 
emitted photons. 

The effects on BR(x2 ~^ Xil) of varying all the masses in the sfermion sector is also 
considered. In particular, the characteristics of the top-squark sector are quite relevant for 
the radiative decay width [^, |16|, jl^. In our previous studies 0, 0, we assumed all the left 



and right squarks degenerate in mass, in order to simplify the multiparameter dependence 
of the neutralino phenomenology. We also neglected the effects of a possible top-squark 
mixing. Here, we will examine the general case and we will see that the behavior of the 
BR(x2 Xil) can be affected by the top-squark sector parameters in different ways when 
the radiative BR is enhanced dynamically or kinematically. 

The outline of the paper is the following. In Sec. 2, we review the theory of the radia- 
tive neutralino decay in the MSSM and fix the notations. We also introduce the possible 
scenarios where a large BR(x2 Xil) regime can arise. In Sees. 3 and 4, respectively, we 



^In this paper, we follow the same convention as in Rcf.[|l| for the sign of /i. We also assume Mi^2 > 
and tan/3 > 1. Note that large values of tan/3 ( 60) are disfavored by a radiative electroweak symmetry 
breaking in the MSSM 
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go through the SuSy parameter regions where a dynamical and a kinematical enhancement 
of the radiative decay can take place. In Sec. 5, we perform a numerical analysis of the 
radiative BR in the relevant parameter regions. Finally, in Sec. 6, we study the top-squark 
sector influence on the radiative neutralino decay. In Sec. 7, we draw our conclusions. 

2 Enhanced BR(x2 XiT) regimes 

The radiative neutralino decay receives contributions, in a convenient gauge 0, from 16 
graphs with all the charged (both SuSy and non-SuSY) standard particles flowing in the 
loop. The corresponding Feynman diagrams are displayed in Fig. |l|. As a result, X2 ~^ Xi7 is 




(a)i (c), (e), (g), 
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Gy\ y 



(b)i (d), (f), (h), 

Figure 1: Feynman diagrams for the radiative neutralino decay X2 ~^ Xil^ i^i the gauge of 
Ref. P]. For each graph shown, there is a further one with clockwise circulating particles in 
the loop. 

a very interesting process in itself, since it is influenced, to some extent, by all the parameters 
and sectors of the MSSM. However, being of higher order with respect to the main tree-level 
decays, the channel X2 ~* Xi7 typically has a BR not larger than a few percent. In previous 
studies [jl6|, the possibility that BR(x2 X1I) g^ts large in the special case X2 — 



^5 — 7 has been considered, according to approximate formulas for the matrix element 
of the process, including only the main contributions from the W'^/x'^ and t/ti loops. In 
Ref. after the full calculation of the matrix element and the decay width, two examples 
of scenarios with large BR(x2 X1I) ^''^^ given, in two different limits. Asymptotic values 
of the relevant SuSy parameters, according to which the actual radiative processes are 
7 — > H'j and Hi i^27, are considered (see below). In Ref. 0, by using the full calculation 
and assuming the gaugino-mass unification, we stress the presence of non- asymptotic regions 
of the parameter space, where the radiative process is still enhanced. In the following more 
general approach, we show that all the above scenarios are just particular realizations of the 
two main enhancement mechanisms for the radiative decay BR. 
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The naive expectation that the BR for X2 ~^ Xi7 is neghgible with respect to the BR's for 
the tree-level neutralino decays is not realized whenever the latter channels are suppressed 
for reasons that either do not affect or affect to a minor extent the radiative process. This 
can happen basically in two cases: 

a) Dynamical suppression. 

Neutralinos are in general superpositions of gauginos (7 and Z) and Higgsinos {Ha and 
The couplings of sfermions to neutralinos involve only the gaugino components 
[apart from terms ~ {mf/Myy), where is the mass of the standard fermion also en- 
tering the vertex], while the Z° only couples through Higgsinos jl^. Then, the direct 
tree- level decays X2 ~^ Xiff require either simultaneous gaugino components in both 
Xi and X2 (for the sfermion-exchange process) or simultaneous Higgsino components 
(for the Z°-exchange process). This is partly true even when the exchanged particles 
are on their mass shell, that is when the two-body channels X2 ~^ ff and/or X2 Z^Xi 
are open. The above requirement does not hold for the radiative decay, since in general 
both the gaugino and Higgsino components of neutralinos are involved in each graph 
of Fig. D, apart from the diagrams (a)-(b) for a massless fermion /. Hence, whenever 
Xiixl) is mainly a gaugino while X2iXi) is dominated by the Higgsino components, the 
tree-level X2 "width for direct decays falls down and the BR(x2 ~^ X1I) is enhanced. 
For pure gaugino/Higgsino states, BR(x2 Xil) can reach 100%. In particular re- 
gions of the parameter space, this picture can be modified by the presence of a light 
chargino. Indeed, the cascade decays X2 Xii^ Xi/i/i)/2/2) when kinematically 
allowed, can take place even for different physical composition of the two neutrali- 
nos, through l^*-exchange graphs involving both W^'s and H's. Note that a different 
dynamical suppression (that we call reduced dynamical suppression) of the tree-level 
direct decays can take place when the sfermions are heavy and the Z°-exchange chan- 
nel is dynamically suppressed by the presence of a dominant gaugino component in 
at least one of the two neutralinos. The latter case will be of some relevance in our 
following numerical analysis. 

As for the decays into Higgs bosons, when m/jo(m^o) < {m^o — rn^o), the two-body 
channels X2 ~* Xih°{A°), open up too. Naively, the latter do not seem to suffer from 
any dynamical suppression, since the MSSM predicts the vertex ZHih°{A°). However, 
we will see that an effective dynamical suppression can be achieved when one of the 
two lightest neutralinos is dominated by a 7 component (not just any gaugino!). In 
this case, due to the absence of the '^Hih'\A°) vertex, the neutralino decays into Higgs 
bosons are depleted as well, and the BR(x2 ~^ Xil) can still be non- negligible. 

b) Kinematical suppression. 

When X2 ^^cl Xi tend to be degenerate in mass, the widths for the different X2 decay 
channels approach zero differently as the quantity e = (1 — m^o/m-j^o) vanishes. For 

^For the neutralino/chargino sectors, we use notations similar to Refs.|l^ [l^. In particular, for the 
neutralino mixing matrix, we use the basis (7, Z, Ha, Hi,), instead of {B, W3, Hi, H2), used in Ref.j^]. This 
choice is particularly suitable for our purposes. 
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the radiative decay, one has P 

where gyfiyp^-y oc 6(7^/167?^ is an effective couphng arising from the one-loop diagrams in 
Fig. 1^ (in general a complicated function of all the masses and couplings to neutralinos 
of the particles circulating in the loops). 

On the other hand, the three-body direct tree-level decays receive contributions from 
either Z^-exchange graphs or sfermion-exchange graphs. The former ones, involving 
the Higgsino components only, in the limit of small e and massless standard fermions 
/, lead to a total width of P, IT] 



Y.nZ - X°l//)zO-e.ch. ^ , (2) 

J ^ ^■O TT IVl ^ 

where Cw ~ 10~^ and C/^ is a number < 1, depending on the Higgsino content of 
the neutralinos (for pure Higgsinos: ffj = 1). Equation (|^) implies a sum over colors 
and five (six) flavors of final-state quarks (leptons). A similar behavior is found for 
the sfermion-exchange graphs.[] One has, for a single channel into a given // pair, 
mediated by a left- or right-sfermion 



r(x°2 - fi//)/_c.. (3) 

fL,R 



ml 



where C^'^ is typically ~ 10^^ — 10^^, but can be slightly larger or much smaller, 
depending on the gaugino content of the neutralinos and on the specific channel con- 
sidered. Yukawa couplings of the Higgsino components to //, as well as L-R mixings 
for the exchanged sfermions /, are here coherently neglected, since we work in the 
massless / approximation and top (s) quarks are not involved in the problem, rrit being 
too heavy. As for the interference term of the Z°- and sfermion-exchange graphs, we 
expect of course the same fifth-power behavior as in Eqs.(|^) and (|^). This implies 
that the ratio of the direct tree-level and the radiative decay widths tends to vanish as 
e^, when e — > or m^o — > m^o. Note also that the presence of the small number C^, 
in Eq.(|D can partly compensate the additional factor of order agm in Eq.(p. Hence, 
for neutralino masses of interest for present/near-future collider physics and for this 
work (where rn^^o^ ~ Mz), when {m^o — m^o) ~ 10 GeV it is already possible for the 
radiative decay BR to receive a substantial factor of enhancement ~ 10^, especially if 
the sfermion-exchange channels are suppressed for some dynamical reason. 

At this point, it is important to stress that, when neutralinos are degenerate within less 
than about 10 GeV, the asymptotic formulas given above for T{x2 ~^ Xif f)z"- /-exch 



^ As can be inferred, for instance, from the treatment of analogous channels for the gluino decay to a 
photino plus a quark pair, in Ref. 0|. 
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are no longer a valid approximation when ruf cannot be neglected with respect to 
{m^o — m^o), as in the case, e.g., f = b or T.f\ 

However, in the subtle kinematical regions around the various heavier rrif "thresholds", 
the X2 ~^ Xiff channels with lighter /'s will dominate over the ones with f = b,T, . . .. 
Hence, neglecting rUf everywhere and using Eq.(^) and Eq.(^) summed over all non- 
top flavors, will have the net result of an overestimate of the total width for the 
X2 tree-level decays and an underestimate of BR(x2 ~* Xil)- The latter simplified 
treatment of the problem (which is the one we will adopt in all the following numerical 
analyses), however, is here justified by at least two good motivations. First, we have 
already seen that for rn^o^ ^ Mz and when m^o — m^o ~ 10 GeV ~ 2m6, that is 
close to the heaviest ruf "threshold", the radiative decay BR can already receive a 
^ 10^ enhancement factor. This factor can be extracted directly from Eqs.(0), (^, 
which in this region are still a reliable approximation. Second, considering values for 
neutralino masses and the difference {ttl^o — m^o) with a precision at the level of a few 
GeV or less, for a given set of input parameters, is not fully sensible when neglecting 
radiative-correction effects on the spectrum, as we do in this paper. (A more extensive 
discussion of this problem can be found at the end of this section.) 

Having the above caveats in mind, we will not be concerned in the rest of the paper 
with fine behavior and subtle kinematical effects for neutralino mass differences at the 
level of several GeV or less. Nevertheless, we expect that the numerical analyses to 
be performed in the following are valid in most of the interesting cases. Even for very 
small — m^o), our identification of the parameter space regions where this can 

happen, and the result of a large BR(x2 ^ Xil) (with approximately the numerical 
value we will indicate) keep holding. 

Regarding the cascade decays through hght charginos, they are at least as kinematically 
suppressed as the normal direct three-body decays. Indeed, when e ^ 0, the width of 
each of the two steps of the cascade decay has a fifth-power asymptotic behavior for 

ei,2 0, where ei^2 = (1 ~ ^x^/^x'i)^ ~ '"^xll'^^x^^- Furthermore, ei^2 < e, since 
m^o < m,-± < m^o must hold, for the cascade to take place. On the other hand, some of 

Xi ^2 

the channels for this class of decays will not be suppressed by possibly small couplings 
like the C"s, in Eqs.(^ and (0). As for situations with very small neutralino/chargino 
mass differences, similar remarks as for the direct case above hold here. 

Finally, the two-body decay into Higgs bosons cannot take place when the mass dif- 
ference between the two lightest neutralinos is less than a few tens of GeV, because of 
the current experimental limits on m^o and m^o 



The conditions a) and b) can be translated into requirements on the SuSy parameters 



^One should, instead, use the full formulas, including finite to/ effects, to get the correct trend of the 
width as a function of the mass difference Am = m^o — to^o — 2to/ = to^o(1 — e/) between the initial 
and final states, as it approaches (or gets even smaller than) m/. One expects that the simple fifth-power 
behavior in terms of e ~ e/ in Eqs.(||) and (^ will be gradually spoiled while entering the "intermediate" 
kinematical region where e/ becomes sensitively different from e. Here, one would have to consider much 
more involved formulas for the widths as functions of e/ . Eventually, when Ato becomes much smaller than 
mf too, a simple asymptotic regime is again effective, but in terms of ey, and possibly different from the 
former fifth-power-like one. 
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tan/?, /i, Ml and M2, which set the mass matrix of the neutrahno sector. 
The tree-level neutralino mass matrix reads, in the convenient basis 
[— ^7, —iZ, Ha = Hi cos (3 — H2 sin (3, = Hi sin [3 + H2 cos /3] 



Ml cos^Ow + M2 sm^Ow (M2 - Ml) sin 9^ cos 6*^ 

(M2 — Ml) sin 9w cos 9w Ml sm^9^- + M2 cos 9y^,- Mz 

Mz sin 2/5 -/i cos 2/5 

-/X cos 2/5 -/i sin 2/3 



(4) 

It is easy to recognize in Eq.(^ two 2x2 blocks, which correspond to: 

(i) the gaugino mass terms, parameterized by Mi and M2 and mixed by the weak angle. 
Their source is soft SuSy breaking. 

(ii) the Higgsino mass terms, parameterized by yU and tan/5 = ^2/^1, whose source is a SuSy 
term in the MSSM Lagrangian, which mixes the Higgs doublets. 

Then, there are only two off-diagonal entries non included in the two 2x2 blocks, corre- 
sponding to Ha-Z mixing terms and equal to M^, which come from the HHZ couplings and 
the SuSy Higgs mechanism. As a consequence, apart from the asymptotic cases where Mi^2 
and/or are much larger than the Z^ mass, it is not possible to have either a pure Z or a 
pure Ha- Hence, whenever a neutralino has a sizeable Z component, it must have a sizeable 
Ha component as well (and vice versa). This means that a neutralino can be a pure gaugino 
only when it is a photino, and a pure Higgsino can only be of the Hi, type (sometimes called 
"symmetric Higgsino", with notation Hs)- Note also that when Mi = M2 (or tan/5 = 1) 
the off-diagonal terms within the 2x2 gaugino (or Higgsino) block disappear. The limits 
(Ml — M2) and tan /5 — > 1 will be crucial for the enhancement of the neutralino radiative 
decay. 

The outcome of the neutralino mass matrix (El) in terms of the neutralino physical com- 



positions and mass eigenstates has been extensively studied in Ref.|jT8|. We use the results of 
that analysis and concentrate here on what is relevant in order to realize either a dynamical 
or a kinematical enhancement of the X2 ~^ Xi7 decay. 

As already mentioned, some asymptotic regimes, where X2 ~^ Xi7 is enhanced, were 
anticipated in Ref.|^. Some enhancement is expected in the following two cases: 

(i) Light-Neutralino Radiative Decay (Mi,M2, < Mz). 

Then, a dynamical X2 ~^ Xi7 enhancement is realized, since Xi 7 X2 ~^ Hb or 
vice versa. In Ref. [0, the limit = 0, 7 — > H'y is treated analytically. Such small 
values of the parameters are not yet excluded by LEP data, provided 1 < tan/5 ^ 2 
(what is called sometimes "light Higgsino-gaugino window" |^). In fact, in this par- 
ticular region, a number of things happens: 

(a) when \fi\/Mz, M2/MZ 0, the chargino mass generally satisfies the current LEP 
lower bounds; 

(b) if also Mi/Mz — >■ 0, the neutralino mass eigenstates are: 7, Hb (with mass eigen- 
value — > 0) and the symmetric and antisymmetric combinations of Z and Ha (with 
mass —>■ Mz). Then, when tan/5 — > 1, the light Hb decouples and the neutralinos 
can only interact with the Z° boson through the vertex Z^HaHb, which is largely sup- 
pressed by the phase space at LEPl energies, since {mj^^ + ruj^J ^ Mz- 

As a result, the data on the Z° peak can hardly constrain this particular region. 
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For instance, in Ref.[^, the analysis is performed by considering the bounds on 
r{Z° — >■ xlxi) from the invisible Z° width and those on r(Z° — XiXh X2X2) froin 
the direct search of neutralinos. However, only the decays X2 ~^ Xi^^*^ ~^ Xiff 
fully taken into account, while the radiative decay, although generally dominant in 
this region, is not properly stressed. In the analysis of Ref.||20||, a tighter BR bound 
from the radiative neutralino decay has been included as well as the effects of data 
taken at ^/s above Mz- In spite of that, part of the region in the {fi, M2) plane with 
M2, ^ 10 GeV still survives for < and tan/3 close to 1. 

This region is wider when Mi is taken as a free parameter and allowed to be quite larger 
than M2. Some significant improvements in probing the above region could come from 
a careful analysis of the data of the recent LEP short runs at ^/s = 130-136, 161 and 
172 GeV and the future runs around ^/s = 190 GeV. If the light gaugino-Higgsino 
scenario is realized, all charginos and neutralinos are expected to have masses in the 
kinematical reach of LEP2 and should not escape detection. 



(ii) Higgsino-to-Higgsino Radiative Decay (\n\,Mz <C Mi,M2 ~ TeV). 

This corresponds to a particular asymptotic case of the kinematical X2 Xil en- 
hancement. Indeed, in this situation the two lightest neutralinos have nearly degen- 
erate masses close to and are both almost pure Higgsinos. Hence, the direct 
tree-level decays can only proceed through Z^-exchange graphs and the ratio between 
the corresponding width and the Hi 2 -f^2,i7 width can be obtained from Eqs.(0) 
and (|1|), and is independent of the sfermion masses. In a sense, this is an optimization 
of the kinematical-suppression mechanism, since such a dynamical suppression of the 
contributions from the sfermion-exchange diagrams allows the fifth-power asymptotic 
behavior of the tree-level decay width and the enhancement of BR(x2 Xil) to get 
always effective for values of {fn^o — fn^o) at the level of 10 GeV or less [cf. Eqs.(l-3)]. 
On the other hand, in the present case, the factor in the numerator of Eq.(|^) 
is close to its maximum 1, hence depleting a bit the radiative-BR enhancement [cf. 
Eq.(H)]. In the asymptotic limit M, one finds 



T{Hi ^ H2ff) 



M1M2 



1 2 



/i(Mi + M2 tan^Oy 



(5) 



where C is a number of order unity weakly dependent on the ratio M^/yU^. Here, 
for large Mi^2, values of ^ Mz/2 are generally already excluded by LEP1-LEP1.5 
data, since they lead to chargino masses lighter than 45-50 GeV. Hence, Eq.(^ tells 
us that one needs very large values of Mi 2 to get a significant enhancement of the 
radiative decay through this mechanism. In addition, in this region, X2 '^^^ often 
decay through cascade channels into a lighter chargino with a non-negligible branching 
fraction. Numerically, by using the full formulas, we checked that BR(x2 Xil) 
^ 10% always for Mi^2 ^2 — 3 TeV, and ^ 45-50 GeV. Furthermore, it turns out 
that to have a sizeable radiative decay in this case (that is to suppress both direct and 
cascade tree-level decays) , one always has to enforce the condition m^o — m^o ^2 — 3 
GeV, which drastically restricts the photon energy and also corresponds to a critical 
kinematical region, as discussed above. 
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Also note that, for nearly degenerate Higgsinos, the radiative corrections may actually 
spoil the enhancement mechanism or, at least, render the tree-level analysis rather 
inaccurate, even if one takes into account finite rrif effects. For instance, the Higgsino- 
mass splittings in this region can receive radiative corrections as large as their tree-level 
values [i.e., up to ±(5-10) GeV], if the mixing between the two top-squarks is large 

Now, we want to extend the quoted studies by analyzing the more general framework 
where either a dynamical or a kinematical X2 X1I enhancement can be realized, without 
assuming any particular hierarchy between the SuSy parameters and relaxing the usual 
unification condition on the gaugino masses. We will not concentrate on particular limits 
of the SuSy parameters such as Mi and/or M2 — >■ and/or — > 0, since they have been 
either already excluded by LEP data or discussed above. 

In the following, we will neglect the effects of the radiative corrections on the neutralino 



mass matrix elements. In Ref . ||23|| , the full calculation has been carried out. Regarding 
the radiative corrections to the neutralino mass eigenvalues, they are found to be generally 
at the level of 3-8% and of the same sign (positive) for all the mass eigenvalues. Only 
occasionally, the lightest neutralino mass can receive larger corrections. No conclusion can 
be easily extracted from that analysis about whether or not and how much the radiative 
corrections may change the composition of a neutralino eigenstate and, in particular, to 
what extent, for instance, an almost pure photino at the tree level, could turn out to be 
a more mixed state at the one-loop level. Our scenarios for an enhanced radiative-decay 
regime rely on a certain amount of adjustment between different SuSy parameters in order 
to get either pure compositions for the eigenvectors or degeneracy for the eigenvalues of 
the neutralino mass matrix. As for the kinematical-enhancement mechanism, we will show 
that a mass difference m^o — m^o ~ 10 GeV, with fn^o^ ^ Mz, is in general sufficiently 
small to get a BR(x2 Xil) of order 40% or more. Since the higher-order corrections 
to the neutralino masses have generally a fixed sign, one can expect a common shift of 
the different masses, while the relative mass differences change only slightly. Thus, we 
expect that the kinematical suppression keeps almost unchanged for m-j^o — m^o ~ 10 GeV 
and that our treatment substantially holds even after radiative corrections, with a possible 
slight redefinition of the interesting regions in the SuSy parameter space. Concerning the 
dynamical-enhancement mechanism, similar arguments may be used, although in this case 
there are less clues to guess the effects of radiative corrections. However, we will see that 
the amount of parameter adjustment required for the mechanism to be effective is not very 
large. For instance, there are significant regions in the parameter space where |(Xi2l-^6)P 
and/or |(x2il7)P ^^^?/ about 0.8 and the dynamical-suppression is still effective, with 
the BR(x2 Xil) of order 50% or more. Hence, an adjustment of the parameters at the 
level of 20% should survive the inclusion of the radiative corrections. 

An additional remark is due for the case of the kinematical enhancement mechanism. 
Since the latter arises from situations where the two lightest neutralinos are close in mass, 
a too strong degeneracy may prevent the experimental detection of the X2 decay, due to the 
emission of too soft photons. In general, one can ensure the presence of a useful experimental 
signature and the phenomenological relevance of the neutralino radiative decay by requiring 
m^o — m^o ^ 10 GeV. This, of course, effectively depletes the actual BR(x2 Xil) that 
can be achieved by the kinematical mechanism in a real experimental framework. On the 
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other hand, if the available cm. energy is large enough, the X2 can receive a sizeable boost, 
even when {m^o — rriy^o) is very small (for instance, when produced in association with a Xi 
at LEP2, or with a Xi at the TeVatron). In order to assess to what extent this is true, one 
can take into account that, assuming an isotropic radiative decay of the produced xl (that 
is neglecting spin-correlation effects), the resulting photon has a flat energy distribution in 
the laboratory, with end points 



mm, max 



\ 



El 

^1 



m-o — m-o 
2m^o 



where E2 is the production energy of the X2 laboratory. 



3 Dynamical X2 ~^ X1I enhancement 



0. 



As already seen, requiring one pure gaugino and one pure Higgsino eigenstate from the 
matrix (D implies both 

Ml = M2 (6) 



and 



tan/3 = 1 . 



(7) 



One then has a pure 7 with mass Mi(= M2) and a pure Hb of mass —fi in the neutralino 
spectrum. In this limit, the other two neutralinos are mixtures of Ha and Z with mass 
eigenvalues (including their signfl |TH 



m 



(±) 

Ha-Z 



M2 + /i± V(M2-^)2 + 4M| 



and mixing angle 







1 COS0 


} V2 



1 ± 



M2 -/i 



(M2-/i)2 + 4M| 



1/2 



(9) 



Requiring that the pure states correspond to the lightest neutralinos, Xi and X2) the absolute 
values of both the eigenvalues in Eq.(|D have to be larger than both Mi (or M2) and In 
the parameter space not yet excluded by the LEP data, this can be achieved only if 



fi<0 . (10) 

Indeed, when fi is positive, the smallest absolute value in Eq.(H) corresponds to choosing 
the negative sign before the square root. It is then sufficient to look at m^l^. The latter is 
always smaller (greater) than both ^ and M2, whenever ^ and M2 are larger (smaller) than 
Mz/2. On the other hand, if /i < Mz/2 < M2 or M2 < Mz/2 < /i, m^^l^ can still (but 

^Of course, the physical neutraUno masses are always positive, but the sign of the mass eigenvalue has its 
own physical meaning, being connected to the neutralino CP quantum numbers and entering the expressions, 
in the basis we adopt, of the Feynman amplitudes for processes involving neutralinos 
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not necessarily has to) be large enough to allow the mass ordering needed for the dynamical 
suppression. However, values of M2 and/or \fi\ ^ Mz/2, with fi positive, are generally 
excluded by LEP data. This is true either because of the chargino-mass bound, or because 
of the direct searches of neutralinos, even without gaugino-mass unification assumptions, 
and even in that window with very small and M2 we treated above, for any tan/5. As a 
result, the dynamical enhancement can be present only for /i < 0. Note that the eigenvalue 
corresponds to a massless neutralino when /x = M2 = Mz- 

Quite different is the situation for /i < 0. Points in the SuSy parameter space where 
f^^H-z both heavier than and M2 do exist for small values of and/or M2. This is 
also true for large values of and M2, i.e. far away from the LEP exclusion region. For 
instance, let's examine the case /i = —Mz- In this case, whenever < M2 < Mz/2, Im^fflzl 
is always less than |//| and corresponds to the next-to- lightest neutralino, the lightest one 
being the photino with mass Mi = M2. Then, no dynamical enhancement can take place. 
For Mz/2 < M2 < Mz{l + \/3)/2, the neutralino eigenstates corresponding to m^l^ are 
always the two heaviest ones. Then, the radiative decay Hb 77, for M2 < Mz, or 
7 — i> 7-fff,, for M2 > Mz, benefits from a strong dynamical suppression of the tree- level 
decays. When M2 gets larger than 1.37Mz, tti'h-z becomes the next-to-lightest neutralino, 
while Hb is the lightest one. Once more, no dynamical enhancement can occur. This rather 
complex behavior generates sharp and well outlined contours for the BR(x2 x/il) in the 
(Ml, M2) plane, in the vicinity of the diagonal Mi = M2 (see Sec. 5). 

Of course, when (Mi — M2) and (tan/3 — 1) go away from 0, all the arguments given 
above, including the condition (|l^), have to be intended in a weaker sense. Large radiative 
BR's can be obtained, even in points of the parameter space where |Mi — M2I ^ Mz/2 and 
tan [3^2. As for the validity of the condition (|T0|) in a less restricted case, one has to note 
that a massless (or very light) state is present whenever the equation 

MiM2/i = (Ml cos^ 9w + M2 sin^ 9w)Ml sin 2(3 (11) 



approximately holds ||18|. This implies a positive /i. For general values of Mi 2 and tan/3, 
this state is a superposition of all the four interaction eigenstates. When (Mi — M2) 
its photino component tends to vanish, while its component gets smaller and smaller as 
tan /3 — > 1. Hence, in general, in a wide region in the vicinity of the curve defined by Eq.(|TI|) 
in the plane (Mi,M2), the lightest neutralino (or the next-to-lightest one) is a mixed state 
made of Higgsino-Z-ino and no dynamical-enhancement mechanism can be present, even if 
Ml ~ M2 and tan/3 ~ 1. 

In our numerical analysis, we have not found any case of sizeable dynamical suppression 
for positive values of /x, in the allowed regions. Note that, for a positive /i, a very light 
chargino can be present. In particular, if 

M2/i = M^sin2/3 (12) 

the light chargino is massless. This gives rise to "forbidden" regions in the SuSY-parameter 
space in the positive-/! case, where our basic assumption LSP= Xi "^^J iiot be satisfied. 
On the other hand, even if m-± > m^o, the chargino often turns out to be lighter than the 
X2-I opening the cascade channels, which do not suffer from dynamical suppression. This 
happens, for low tan/?, only rarely for < and often for /i > 0, and provides a further 
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explanation for the lack of significant BR(x2 ^ Xi7)"6nliancement regions in the positive-/^ 
case. 

In order to get a general insight of the neutralino physical composition pattern, we show 
in Figs. 2-4 the behavior in the (Mi,M2) plane of the quantities A = (Xil7)^(X2l-^f))^ 
B = {Xi\Hb)'^ {X2\^)'^ (^) -S < 1), which express the neutralino physical "purity" in the 
dynamical enhancement framework. We also study the effect of varying both the fi and 
tan/3 values in the interesting ranges. In the following, we will see that, in order to achieve 
an appreciable dynamical enhancement of BR(x2 Xi7) 5 either A or B should be as high 
as 0.8-0.9. This condition is fulfilled in a substantial portion of the (Mi,M2) plane, for 
low tan/5. As expected, considerably high values^ of A, B can generally be achieved for 
Ml ~ M2, when tan/5 is close to 1. Some deviations from the expected behavior in the 
limit Ml = M2 and tan/? = 1 we treated above are due to our choice tan/3 = 1.2. In 
general, the presence of contour lines that delimit an abrupt change of regime in either A or 
B generally corresponds to a crossing in the mass ordering of a physically "pure" neutralino 
and a "mixed" state (or of two "pure" states). 




Figure 2: Contour plot for the quantities: (a) A = (Xi|7)^(X2l-^b)^ (b) B = 

{xl\Hby{X2\l)^ (AB < 1), in the case tan/? = 1.2, fx = -2Mz. A and B give a 
hint of the "purity" of the limit Xi = "1^ X2 — -Bb or vice versa. The line corresponding to 
gaugino-mass unification is also shown. 

For instance, in Fig. 2, proceeding along the Mi = M2 diagonal from small to large 
Ml 2 values, one can single out four different regimes and this can be explained by the 
discussion above. The behavior of A and B along this diagonal is of particular interest, as 
anticipated. Indeed, the latter is the only region where A, B can substantially exceed the 0.8 

^°Note that in the "democratic" case (Xi It)^ — iXil^b)^ = 1/4; * = 1, 2, one would have A=B=1/16 only. 
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level and the dynamical enhancement mechanism can be fully effective. When Mi^2 ~ 150 
GeV, the lightest neutralino is an almost pure photino with mass close to Mi 2 and the 
next-to-lightest neutralino is a mixed Higgsino-Zino state with mass close to Eq.(^_|_ (the 
-|- refers to the sign considered in the equation). Also, xl ~ Hf, with mass close to — /i, 
and xl is the other mixed state with mass close to the absolute value of Eq.(|^)_. For 150 
GeV ^ Ml 2 ^ — fi = 2Mz the mostly-i^b and the lighter mixed states exchange their 
role, becoming the X2 Xs; respectively. This arrangement is then suitable for a 

dynamical BR(x2 ~^ Xil) enhancement with xl — 1 and X2 — [cf. Fig. 2(a)]. For 
2Mz ^ Ml 2 ^ 200 GeV, the mass ordering of the dominantly 7 and Hb states is exchanged, 
but one still has a scenario with dynamical enhancement [cf. Fig. 2(b)]. Note also that 
the two mixed states exchange their role as well, the one corresponding to the negative 
mass eigenvalue becoming lighter than the other. This double level crossing takes actually 
place in the highly degenerate point Mi = M2 = —/x sin 2/5 ~ 179.4 GeV, where all the 
contour levels in Fig. 2 tend to crowd. Points of this kind will turn out to be of relevance 
for the kinematical mechanism too (cf. Sec. 4). Finally, when Mi 2 ^ 200 GeV the negative 
eigenvalue corresponding to the mixed state crosses the mass level (~ Mi^2) of the almost 
pure photino, which then becomes the Xs- Then, the dynamical mechanism stops working, 
in spite of the presence of an almost pure Hb as Xi- 
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Figure 3: The same as in Fig.0, but for 



-M. 



In Fig. 3, we show how the general picture for A and B evolves when /x goes from —2Mz 
to — M^. One can check that the situation is qualitatively similar to the previous case, once 
the whole structure of the contour plots in the (Mi,M2) plane is shifted toward the new 
crossing point Mi = M2 = — /isin2/? ~ 89.7 GeV. Note, however, that here the quantity 
A can reach higher values (of order 0.9 or more), and the region where B is large is also 
wider in Fig. 3 with respect to Fig. 2. Furthermore, it is interesting that it is possible to 
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approach high values of A and B (A ~ 0.9 and B ^ 0.85 in a certain Mi^2 interval) even in 
the gaugino-mass unification case. 

For tan (3 as high as 4 (Fig. 4), A and B never reach 0.8 and, consequently, never prompt 
a sufficient dynamical BR(x2 ~^ Xil) enhancement. 
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Figure 4: The same as in Fig.H, but for fi = —Mz and tan (3 = 4. 



4 Kinematical xi ~^ Xi7 enhancement 

As anticipated in Sec. 2, when Xi and X2 almost degenerate and [tti^o — m^o) is smaller 
than about 10 GeV, for fn^o^ ^ Mz, the radiative decay is enhanced by a large, purely 
kinematical factor. The latter can actually turn out to be an overall factor, especially when 
the contribution of the sfermion exchange to the X2 tree-level decays is suppressed. Hence, 
one has an optimization of this "kinematical" enhancement for heavy scalar masses and/or 
small gaugino components in Xi and/or X2j as can be inferred by comparing Eqs.(0) and (^. 
Indeed, after summing over ffavors and colors in Eq.@, if (some of) the C^'^ are not far from 
1 and/or (some of) the Jl^fi are not much heavier than about Mz, the sfermion-exchange 
channels and the interferences can dominate in this regime. The net effect of this can be to 
slow down or even prevent the introduction of an effective kinematical enhancement, so that 
BR(x2 ~* X1I) might get large, if ever, only for neutralino mass differences smaller than a 
few GeV, which we know to be a potentially dangerous region to explore, at least in our 
approximat ions . 

Given the above considerations, the next step is to find out where in the SuSv-parameter 
space the masses of the two lightest neutralinos are almost degenerate. To this aim, one has 
to consider the 4*'^-degree eigenvalue equation associated with the neutralino mass matrix (^) 
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= m- + Am^ + Bml + Cmi + D 
with : 

A = -Tr(A^^o) = -(Ml + M2) , 
B = M1M2 - ii^-Ml , 
C = 

D = det{M^o) = fi[Ml{Mi cos^^^ + M2 sin^Ow) sin 2/3 - /iMiMs] , (13) 



IVI1IVI2 — /X — IVl^ , 

fi\Mi + M2) + M2(Mi cos^^^,, + M2 sin^^^^ - /i sin 2/3) , 
det{M^o) = /i[M|(Mi cos^^w + M2 sin^Ow) sin 2/3 - /iMiMs 



where (i = 1, . . .4) is the general neutrahno mass eigenvalue. Then, one has to force 
Eq.(|l3]) (which, arising from an Hermitian, real and symmetric matrix, has four real roots) 
to have (at least) either two identical roots or two opposite roots, m° = ±m° (when using 
the superscript °, we generally refer to a degenerate eigenvalue ).|^ 

Exact expressions for the neutralino masses and mixing can be found in Ref.[^]. Here, 
we are mainly interested in special cases for which it is possible to extract approximate 
formulas, more useful for a physical interpretation. 

Involving a linear combination of 12-dimensional terms (where up to the 4*^^ power of 
one of the coefficients A, B, C, D can appear), the general necessary condition to get two 
identical roots from Eg. ([T3|) , seems far too complicated to give any useful information and 
even to be displayed here. As for the case of two opposite roots, a simple necessary condition 
in terms of A, B, C, D can be derived: 

A'^D + ABC + C'^ = . (14) 

Unfortunately, Eg. (|T^) , when translated in terms of Mi, M2, fi, sin 2/3 and sin^6',^, turns out 
to be guite complex too. Thus, in the following, we consider only interesting limits of the 
SuSy parameters, such as tan/3 1 or Mi —>■ M2. In this way, by reducing to zero the 
off-diagonal terms in one or both of the 2x2 blocks of the matrix (^, one considerably 
simplifies the eigenvalue eguation (0) and allows to disentangle the relevant degeneracy 
scenarios. In this way, we will single out some sufficient conditions, which ensure exact 
or approximate degeneracy between the two lightest neutralino states. This procedure will 
be supported by an extensive numerical study (scanning the whole SuSv-parameter space) 
of the relevant mass splitting {m^o — m^o), which is shown in the following, as well as a 
numerical analysis of BR(x2 Xil)^ that we present in the next section. On the basis of 
the numerical analysis, we found that necessary conditions for an exact degeneracy of the 
two lightest neutralino states are 

tan/3 = l or Mi = M2 , (15) 

with the additional reguirement yU < 0. The latter directly translates in some necessary 
conditions for a sizeable kinematical enhancement of BR(x2 ~^ Xil)- Interestingly enough, 
these conditions are the same as the ones we found in Sec. 3 for a dynamical enhancement, 
although in the latter case, the two reguirements on tan (3 and Mi 2 have to be fulfilled at the 
same time. The conditions ([T5|) for a BR(x2 Xil) kinematical enhancement are of course 

^^The physical neutralino mass is given by m^a = \mi\. 
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interesting only in regions allowed by LEPl-1.5 data. Furthermore, we do not consider here 
asymptotic regimes with M12 ^ Mz and/or ^ Mz, or Mi^2 ^ Mz and/or -C Mz- 
Indeed, for instance, the limit <^ Mi 2, Mz leads to two almost degenerate Higgsinos- 
like neutralinos, as discussed above and in Ref. 0. Other asymptotic cases with possible 
degeneracy have been mentioned in Sec. 2. Here, we will limit our analytical and numerical 
analyses to the region where \fi\ and Mi^2 fall both in the interval [M^/4, 4M^]. 

We stress that here we are interested to the degeneracy of the two lightest neutralinos. 
This singles out only a few among the possible degeneracy scenarios for the four neutralinos, 
and makes the analysis more involved. 

In order to find sufficient conditions for the degeneracy of the two lightest neutralinos, 
let's now first consider the limit Mi = M2. Then, it is convenient to solve the equation ([T3| ) 
just in terms of M2 as a function of the generic eigenvalue rrii. For rrii 7^ ±yU, one then gets 
two branches 



(M2)+ = rrii, 

, ,9"^i + Msin2/3 , , 

M2 _ = rrii -Ml 2 2 • 16 

rrii ~ 

Note that the branch (M2)+ describes the behavior of one neutralino mass eigenvalue only, 
while the branch (M2)_ is threefold and corresponds by itself to three generally different 
eigenvalues. The two main branches (M2)± intersect in a point (corresponding to a degen- 
eracy with same-sign eigenvalues: m^^- = M2), whenever 

Ml = M2 = -/isin2/? . (17) 

This case is quite interesting. Indeed, one has: 

m\ = m°^ = Ml = M2 = -/isin2/5 , (18) 
-m° = m° = vWm|, (19) 

that is, the two mass eigenvalues with the lower absolute values coincide, while the two 
eigenvalues with the higher absolute values are opposite. As for the composition, one of the 
two light degenerate states is a pure photino, while the other is a mixture of Z and Ht, with 
(Zlxg)^ = /i cos ^2,(3/ \J Mz + /i^ cos 2(3. Thus, the scenario of Eq. ([T7|) is relevant both for the 
kinematical and for the reduced dynamical enhancement (cf. Sec. 2), since the Z°-exchange 
contribution to the tree-level decays is highly suppressed. The degeneracy corresponding to 
Eq.(p!8D is of course removed when Mi and M2 get far apart. Nevertheless, one can check 
that for Ml — M2 ~ 10 GeV the degeneracy may still be effective for a sizeable kinematical 
enhancement of the radiative decay. 

In principle, one can get other degeneracy scenarios by using Eqs.(|16]) in different ways, 
but we will show in the following that they are not relevant for our purposes. First of all, it is 
clear that the degeneracy between the heavier neutralinos corresponding to Eq.(|19]) cannot 
be a direct result of the intersection of the two main branches (M2)±. However, it can happen 
that two of the three sub-branches of the branch (M2)_, for a given value of M2, correspond 
to two opposite eigenvalues, and hence to degenerate neutralino masses. This is the case 
when (m^ — yU^) = M|, which is nothing but the degeneracy ([T9|). The additional fact that in 
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such a scenario one necessarily has also M2 = — sin 2f3 [corresponding to a real intersection 
of the two main branches (M2)± and to the further degeneracy (|18D of two always lighter 
neutralino states] makes the above circumstance negligible for us here. One can look for 
other cases of degeneracy by considering the possibility that one of the eigenvalues described 
by the branch {M2)- has the same absolute value of the one corresponding to the branch 
(M2) + , but opposite sign. This corresponds to solving the equation 



2M2 -/isin2/3 



while, for the two degenerate eigenvalues, one has: m° = — = M\ = M2. The solutions 
of Eq. ( pOD are in general complicated expressions, but one can easily find them numerically. 
After such an analysis, we did not encounter any further case of exact neutralino mass de- 
generacy from Eq. ([2D|) in the limit (Mi — M2) — > relevant for the kinematical enhancement, 
in regions allowed by the LEP data. 

Relaxing the limit Mi = M2, a different necessary condition to get mass degeneracy 
in Eq.([T3|) is, indeed, tan (3 —>■ 1. In this limit, contrary to the previous case, the easiest 
part of the neutralino mass matrix is the Higgsino sector. Then, we can solve the eigenvalue 
equation with respect to as a function of mj, in order to get other sufficient conditions, and 
scenarios of interest for the kinematical enhancement. One then finds again two branches 

= -rui , 

2 - Ml cos'^e^ - M2 siir^e^ 
/i + = rrii- Mz —- — , 21 

[nii - Mi)(mj - M2) 

where, similarly to the case of Eqs.(p!6|), the branch describes a single neutralino mass 
eigenvalue and the branch (yu)+ corresponds to three different eigenvalues. Note that here, 
contrary (and complementary) to the case of Eq.(|l6|), one can have rrii = ±/x, but not rrii = 
Ml or M2. By using Eqs.(^ID, one can single out degeneracy scenarios with m° = ±m° = /i. 
In order to realize the case with two same-sign degenerate eigenvalues, the general condition 

/. + Micos^^ +M2sinX_ 
^ (/i + Mi)(/i + M5^ ^ ' ^ ^ 



must hold. Again, the corresponding explicit solutions are rather complex, but one can solve 
Eq.(p2D numerically to find the regions of the parameter space interesting for the kinematical 
enhancement. Since, in contrast to the case of Eq.(^), Eq.(^) gives rise to interesting 
scenarios, it is useful to consider here, in addition to tan/? — >• 1, the special limit siia^Ow = 0, 
which allows a simplified analytical treatment. Indeed, in this way, two free gaugino mass 
parameters are still present, but the mixing in the gaugino sector disappears. Furthermore, 
this limit is not too far away from the real physical case sin^^^v — 0.23. Equation (p2D gives 
then the solutions (when tan/5 = 1) 



= m° = -fi = 1{M2 T \/ + 2M|) . (23) 



2 



The existence of the solutions ( p3l) (not their exact form), is independent of the limit 
siia^Ow — ^ and the corresponding exact degeneracy is removed only for tan/3 7^ 1. The 



S. Ambrosanio and B. Mele / Supersymmetric Scenarios. 



19 



introduction of this limit allows us to disentangle other two interesting cases, as long as 
tan (3 = 1, that is 

m° = m] = -/i = Ml , (24) 

and, when the condition 

7Vf2 

Ml ^ u = (25) 

^ M1-M2 ^ ' 

is fulfilled, 

M^ 

m° = m° = Mi = a ^ — . (26) 

^ ^ M1-M2 ^ ' 

The latter cases, needing mj = Mi, cannot be directly derived from the two branches in 
Eg. (pri) ■ The correct procedure to get them is to solve the eigenvalue equation with respect 
to yU or Ml, by applying both the limits tan/5 1 and sin^6'^ — > simultaneously. 
In order to understand the nature of these additional solutions and their link with the 
limit sin^6'v7 — > 0, some further explanation is needed. A solution corresponding to Eq.(p^) 
survives when sin^^^^ 7^ 0, although the expression for the degenerate mass eigenvalues 
receives some corrections, as in the case of solutions (p^). What makes case ( p^ different 
from the previous ones is that the corresponding degeneracy is not removed when tan /3 goes 
away from 1 and both ^\Y?dyv 7^ and tan /5 7^ 1 are needed to do the job. In this sense, 
the degeneracy corresponding to Eq.(p^ is more "solid" than the others. As for solution 
(p6|), instead, it represents a spurious case which does not correspond anymore to an exact 
degeneracy when sin^^,^ 7^ 0, but only to a case where two neutralino mass eigenvalues are 
close to each other (in the limit tan/3 — > 1), although not quite equal. To get an effective 
kinematical BR(x2 ^ Xi7) enhancement, we are interested in scenarios where neutralino 
mass differences of order 10 GeV or less arise. In order to obtain such a small (m° — m°) 



when Eq.(p6D holds, one needs Mi, M2 -C M^, with ^ TeV [because of condition ( pSf ) 
Then, the quantity (Mi 0.0^0^ + M2 ^Yc?d^) in Eq.(|l3|) can be treated in the same way as it 
would be in the limit sin^^vK 0. Therefore, we will neglect this possibility (corresponding 
to an asymptotic case, already excluded by LEP) and we will focus on solutions (|23|) , (|2^) . 

We stress that the existence of the solutions (|23|) , (|2^) (not the exact expression for the 
degenerate eigenvalues) is independent of the limit sin^6'vi/ — >■ and the corresponding exact 
degeneracy is removed only for tan/? 7^ 1. Anyway, the simplified solutions we found allow 
us to emphasize some remarkable properties which remain valid with a good approximation 
for sin^^vv — 0.23 (and, often, even for tan/? ~ 1, rather than exactly 1). For instance, 
the solutions corresponding to (|23|)+ and ( P^ are only possible for negative values of /i 
[even /i < — M^/y^ in the case (p3D+)].H In contrast, the solution (p3D_ is allowed only 
for < < Mzl\pl (that is in a region that, particularly for small tan/3, is excluded by 
LEPl-1.5 data, due to the presence of a light chargino). Hence, solution (^)_ will not play 
an important role in the following discussions. Also, note that the solution (^^ is present 
irrespective of the particular value of M2 and, similarly, the solutions (^) do not depend 
on Ml, in the limit tan/3 — > 1, ^Yc?dw and exact degeneracy. For ^\Y?dyv ~ 0.23, the 
solution ( p3D [(p^j develops a weak dependence on Mi [M2], as will be shown in the following 
numerical study. 



^The subscript ± of the equation number picks out one of the two possible signs in Eq.d2 
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As for the limit tan (5 1, up to now we only took care of deriving some sufficient 
conditions for the degeneracy of any pair of neutralinos. Now, we need to check when 
the degeneracy scenarios we singled out actually concern the two lightest neutralino mass 



eigenstates. We will focus on the more interesting scenarios (^31)+ and (p^). Which pair of 
mass eigenstates is involved in the degeneracy depends also on the parameters not directly 
entering the approximate conditions (pSf), (p4D, in a generally simple way. For a given value 
of typically one observes that, for Mi [M2] > the degeneracy of the kind (|2^)+ 
[(p4D] indeed concerns the two lightest neutralino states. On the other hand, as long as 
Ml [M2] < the solution for the mass degeneracy corresponds to m^o = m^o = — /x and, 
thus, does not give rise to any kinematical suppression. 

All the tan/3 = 1 scenarios above are derived by forcing the two branches in Eqs.(pID to 
meet in a point of the [rrii, /i) plane, corresponding to a strict degeneracy m° = m°. However, 



as done above in the limit (Mi — M2) with Eg. (pOD , here too we must take into account 
the additional possibility of a degeneracy with opposite eigenvalues m° = —m°. Imposing 
the corresponding condition (/i)+ = — (/u)_ in Eqs. (pi]) leads to the interesting case 

ml = -m] = ^ = Mi cos^e^ + M2 sin^^^^ , (27) 

which satisfies the Eq.(|T^) and is only realized for positive values of /i. This gives rise 
to a scenario where a pure Hi, is degenerate with a superposition of the other interaction 
eigenstates. The other two neutralinos are mixed j-Z-Ha states too, and correspond to the 
mass eigenvalues 

1 



± 

J-Z-Ha 2 



Ml + M2 ± V(Mi - M2)2 + 4M2 



(28) 



On the other hand, the degeneracy must involve the two lightest states. It is easy to show 
with analytical arguments that this can never happen in regions allowed by LEPl-1.5 data. 
Hence, the scenario of Eq . (|28|) is not relevant for the radiative neutralino decay. 

Finally, it is interesting to note that in the limit where both tan /3 — >■ 1 and (Mi— M2) —>■ 0, 
both the kinematical and the dynamical enhancements can be optimized at the same time 
in the special point Mi(= M2) = — /i, where the two lightest neutralinos are always a pure 
7 and a pure Hh [cf. Sec. 3, Eqs.(|D and Eqs.(|T8D, (|19[)]. We will see, in the numerical 
analysis of Sec. 5, that a considerably wide region of the SuSv-parameter space where high 
BR(x2 Xil) values are realized is centered on this highly degenerate point. 

As for the kinematical enhancement, by using numerical methods, we did not find any 
other clear case of X1-X2 exact mass degeneracy, besides the ones we have described above. 
Also, it was not possible to achieve an approximate degeneracy (at the level of a 10 GeV 
mass difference) either, in regions of the SuSy space where the necessary conditions ( [T3| ) are 
rather far from being valid. This does not mean we listed all the possibilities for neutralino 
mass degeneracy. For instance, one can consider the case D = det(7Mj^o) = 0. This can be 
achieved either when fi = and/or Mi_2 = 0, or whenever Eq.(|ll|) holds. The first option 
was already considered among the asymptotic cases. In the second case, one is left with a 
simplified eigenvalue equation, which gives rise to other degeneracy scenarios.0 However, 

For instance, one finds a non-trivial degeneracy, for a given value of tan/3, in the case Mi = M2 — 

f s 1/2 / sl/2 

(Mz/2) n + Vl + 8sin^2/3j and ^ = (Mz/V2) ( VI + 8sin^ 2/3 - Ij . The corresponding degen- 
erate neutralinos are X2 ^-^d Xsi with opposite mass eigenvalues. 
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these scenarios always give rise also to at least a null mass eigenvalue; hence, the degeneracy 
can only concern the X2 ^^e X3 or the X3 the xl- Other complex degeneracy 
situations, not of interest here, can be constructed. 

In summary, the relevant approximate scenarios for the BR(x2 ^ Xil) kinematical 



enhancement are given, for /i < 0, by Eq. ([l8|) , when Mi ~ M2, and by Eqs.(^^+ and 



when tan/3 ~ 1. In order to get a clear picture of the non-trivial behavior of the neutralino 
mass degeneracy, we now show a set of contour plots for the two lightest neutralino mass 
difference {m^o — m.^o) in the (Mi, M2) plane, for different values of fi and tan/3. 




Figure 5: Contour plot for the difference (in GeV) between the two lightest neutralino masses 
in the (Mi,M2) plane for fi = —2Mz, and tan/3 = 1.2 (a) and 4 (b). Different levels are 
represented by lines of different style. The straight line corresponding to the gaugino mass 
unification is also shown. 

In Fig.p, the case /i = — 2M^ is shown for tan/? = 1.2 (a) and 4 (b). The line correspond- 
ing to the gaugino mass unification is also plotted. The general pattern is highly non-trivial 
and quite dependent on tan/3. In Fig.^(a), one can easily note the presence of two quite 
narrow bands, a vertical one and a horizontal one. The horizontal one corresponds to the 
region where the degeneracy scenario of Eq.(PB|)+ is approximatively realized for tan/3 close 
to 1 and is well outlined by the 15-GeV mass- difference contour. As anticipated, the band 
is there only for Mi ^ — fi = 2Mz and shows only a weak dependence on Mi. Indeed, 
by solving numerically Eq. (^2]) for the general condition of degeneracy, we found that the 
contour line of exact m° = degeneracy, in the limit tan/3 = 1, passes through the points 
(Mi,M2) = (200, 168.9); (250, 166.2); (300, 165.7) GeV, showing a small dependence on Mi 
due to the finite value of sin^6',^. Note that the same equation indicates the presence of a 
degeneracy also for M2 = 163.7 and 161.5 GeV, respectively for Mi = 100 and 150 GeV, 
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but this corresponds to a case with 1712 = mg, which of course does not show on the figure 
and is not interesting here. The vertical band, well outlined by the 8 and 15 GeV contours, 
represents the degeneracy (123). Again, the band is present only for M2 ^ '2Mz and has 



only a weak dependence on M2, depending on the finite value of sin^Oyy. The contour line 
of exact degeneracy passes now through the points (Mi,M2) = (179.7,200); (178.2, 250); 
(177.7, 300) GeV. Again, there is also a degeneracy concerning X2 X3 when Mi = 175.7 
and 170.9 GeV, respectively for M2 = 100 and 150 GeV. 

Two additional remarks are in order. First, the vertical band of scenario ( p^ ) is clearer and 
corresponds to an higher degree of degeneracy. This was explained above, in connection with 
the relation of this scenario with the limit siia^Ow 0. Second, both the 8 GeV and the 15 
GeV contour have a bulge where they change direction, along the diagonal, towards lower 
values of Mi and M2. Even more visible is this effect if one observes the 25 GeV contour. 
This corresponds to the region around the point in the (Mi,M2) plane where the further 
case of m° = rrig degeneracy we singled out is realized. This is given by Eq.(|TB|), that, in 
the special case tan/3 = 1.2, gives Mi ~ M2 ~ —/i sin 2/5 = 179.4. Some degeneracy is still 
present along the diagonal Mi = M2, to the left and below the highly degenerate region 
Ml ~ M2 ~ —fi, although to a minor degree than in the horizontal and vertical bands. The 
degeneracy scenario (|18D is more crucial in Fig. ^(b), where tan/5 is far away from 1 and the 
other degeneracy scenarios cannot be realized. Here the 8, 15 and 25 GeV contours surround 
the degeneracy point Mi = M2 = — /isin2/? = 85.8 GeV and extend along the Mi = M2 
diagonal in both the directions. Furthermore, note that this point does corresponds to a 
case of exact degeneracy, while the "median" lines of the horizontal and vertical bands in 
Fig. ^(a) do not, since tan/5 is only approximately equal to 1. It is useful to check that 
the vertical line of approximate degeneracy [see, for instance, the one for Mi ~ 175-180 
GeV in Fig. 5(a)] is directly related to the mass-level crossing of two slightly mixed lightest 
neutralinos along the same line (see, for instance, the behavior of A and B in Fig. 2). 
We will see in the next section that the degeneracy along the Mi = M2 diagonal will trans- 
late in explicit effects on BR(x2 X1I) fo^' tan/? well above 1, while, in the low tan/5 case, 
they will be mixed with and partly hidden by the dynamical enhancement. A final comment 
on Fig. ^ is about the gaugino mass unification. We can see that is not possible in the 
unified case to realize a Xi^X2 approximate degeneracy at a level of less than 25 GeV mass 
difference, unless M2 ^ 300 GeV and tan/? is small (besides the very low Mi^2 region). 

Fig. 6 shows how the general picture evolves when /z goes from —2Mz to —Mz- The behavior 
of {fri^o — m^o) is qualitatively similar to the previous case, although the regions of strong 
degeneracy are shifted towards lower values of Mi^2- Also, it is now possible to achieve 
interesting degeneracy scenarios for which the gaugino-mass unification holds and M1.2 are 
in a range of interest for present collider physics. 

Before coming to the numerical study of BR(x2 Xil)y it is useful to note that the 
dynamical and kinematical mechanism can be present at the same time in special cases and 
work together to enhance BR(x2 ^ Xil)- For instance, a moderate degeneracy is needed 
between m^o and m^o for the kinematical mechanism to be effective, when Xi aiid X2 ^i^o 
have a definite different composition and vice versa. Also note that the general necessary 
conditions (i.e. tan/5 ~ 1, Mi ~ M2, with fi < 0) are the same for both mechanisms, and 
this strengthens the enhancement effect. 
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Figure 6: The same as in Fig. but for fi = 

Numerical analysis of BR(x2 Xil) 



On the basis of the results presented in Sees. 3 and 4 on the enhancement regimes for the 
radiative X2 decay, we are now ready to explain the nontrivial behavior of the corresponding 
BR in the SuSy parameter space. One of the main findings will be the existence of significant 
regions of this space, of interest for collider physics, where X2 ~^ Xil is the dominant X2 
decay. Following the previous discussion, we present the BR(x2 Xil) i^i the (Mi,M2) 
plane, for fixed values of tan/3 and /x. We also discuss the BR dependence on the scalar 
masses. 

First of all, one has to set the parameter regions already excluded by the experimental 
search. We recall that the usual analysis implies the condition M1/M2 = (5/3) tan^ 6'^. 
Relaxing the latter, the definition of the exclusion regions gets of course more involved. 
We considered the following bounds from LEPl data on the Z° line shape and on neutralino 
direct searches [|^, 



BR{Z° 
BR(Z" 



SuSy) < 23 MeV 

SuSy) < 5.7 MeV 

< 3.9 X 10" 

< 3.9 X lO-*^ 



X1X2) 

xlxl) 



6 



where we took into account not only the XiXi contribution to the Z° invisible width, but 
all the channels Z° XaX]^ "^j J = 1; • • • ? 4, with following invisible decays of the produced 
heavier neutralinos Xi ^ Xi^^i- ■ ■ ^^)- 
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As for LEP1.5, we imposed the following limits from the direct searches of neutrali- 
nos/charginos during the runs at ^/s = 130.3 and 136.3 GeV [2G 



< 2 pb at = 130.3 GeV 



if (m^o — m^o) or (m-± — m^o) > 10 GeV 



"X2 



0"v 



< 2.4 pb at = 136.3 GeV , 

if {m^o — m^o) or (^«-± — m^o) > 10 GeV 



which corresponds to not allowing more than ten total visible events from neutralino/chargino 
production at LEP1.5. Here too, the branching fractions into visible final states of the pro- 
duced particles have been taken into account while computing civis- The above LEP1.5 limits 
translate into an approximate bound on the chargino mass m-± ^ 65 GeV, when the sneu- 
trino is not too light (m^; ^ 200 GeV) and there is enough phase space available to ensure the 
presence of rather energetic particles among the chargino decay products: {'>tl^± — m^o) > 10 
GeV. 

A few comments are in order. By applying the above experimental constraints, one finds 
that for rather small values of tan/? and ^ fi ^ Mz wide regions in the plane (Mi,M2) 
are excluded. For instance, at tan/3 = 1.2 and fi = Mz, the area with Mi,M2 ^ 180 GeV 
is forbidden. The regions with positive n are, however, not much relevant for the radiative 
X2 decay, as we know from Sees. 3 and 4. For low tan/5, negative fi and small |//| (around 
50 GeV) one gets highly non-trivial exclusion regions (see below). On the other hand, 
whenever the chargino search turns out to be the most effective means of constraining the 
plane (Mi, M2), the mass limit and the forbidden region tend to get independent of Mi. The 
possible relevance of neutralino searches can translate, instead, in more involved bounds, 
depending on Mi. 

In Fig. 7, we show the BR(x2 Xil) for tan/5 = 1.2 and fi = —Mz, i.e., in a regime 
where both the dynamical and kinematical enhancement can be realized. The shaded region 
is excluded by LEPl-1.5 data. In this figure and in all the following ones, we calculate the 
Higgs sector masses and couplings by assuming rriA = 300 GeV. In Fig. 7(a), the kinematical 
enhancement is optimized by the large value (1 TeV) assumed for the sfermion masses. We 
find a significant area in the plane (Mi, M2) where BR(x2 Xil)> 90%. Its shape can be 
straightforwardly explained by putting together the information from Fig. 3, on the physical 
purity of the two lightest neutralinos, and Fig. 6(a), on the neutralino mass degeneracy. It is 
interesting to note in Fig. 7(a) how much, along the Mi ~ M2 diagonal, different enhance- 
ment mechanisms can be effective and can contribute to a large BR(x2 Xil) ~ 90%, 
depending on the Mi 2 values. For instance, when Mi ~ M2 ~ 40 GeV, the two lightest 
neutralinos have widely different masses [cf. Fig 6(a)] and the kinematical enhancement 
is not operative. Furthermore, since tan/? is not exactly 1, X2 comparable Z and Ht 
components, spoiling the full dynamical enhancement. However, since the sfermions are very 
heavy, it is sufficient the presence of a nearly pure Xi — 7 to largely deplete the only tree- 
level decay channel left (through Z° exchange) and to give rise to very large BR(x2 ^ Xil) 
values {reduced dynamical suppression, cf. Sec. 2). Of course, after lowering the sfermion 
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Figure 7: Contour plot for the branching ratio of the neutrahno radiative decay. In (a), 
all the sfermions are taken degenerate with mass equal to 1 TeV. In (b), the left and right 
charged slepton masses are taken degenerate and equal to 120 GeV; the sneutrino mass is 
calculated by using the SU(2) sum rule; the squark masses are taken all degenerate and 
equal to 300 GeV. The Higgs sector masses and couplings are set by rriA = 300 GeV. The 
shaded region is excluded by LEPl and LEP1.5 data. 

masses, this enhancement regime does not survive [cf. Fig. 7(b)]. Proceeding towards higher 
Mi^2 — —fJ' sin 2(3 values on the diagonal, one has two enhancement mechanisms getting effec- 
tive at the same time. First, the Hb component of the X2 grows [cf. Fig 3(a)], giving rise to a 
"full" dynamical enhancement. Second, m^o and m^o get closer and closer [cf. Fig. 6(a)] and 
the kinematical mechanism becomes effective as well. The most favorable situation is then 
realized in this region. For larger values of Mi^2 in the [100, 120] GeV range, the dynamical 
enhancement dominates, but with Xi — Hf, and X2 — 1 [cf- Fig- 3(b)]. Finally, when Mi 
and M2 are both ^ 120 GeV, large BR(x2 Xil) values cannot be achieved anymore. 

The relative importance of the various kinematical enhancement scenarios can be guessed 
by comparing Fig. 7(a) with Fig. 7(b). In the latter, lighter scalars, that tend to reduce the 
effect of the kinematical suppression in the X2 tree level decays, are assumed. Nevertheless, a 
strong radiative enhancement mainly due to a dynamical suppression of the tree-level decays 
is still present for Mi ~ M2 — —fi. Note that the kinematical enhancement which is expected 
from the scenario [cf. the horizontal band of approximate degeneracy in Fig. 6(a)] is 

only slightly influenced by the change of the sfermion mass values, while the scenario (p^ ) 
(vertical band) corresponds to higher values of BR(x2 Xil) for heavy sfermion masses, 
and is much more sensitive when changing these parameters. This happens in spite of the 
greater degree of neutralino mass degeneracy corresponding to the scenario (p^ . The reason 
is that in the case (|23|), in addition to the kinematical suppression of the tree-level decay 
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Figure 8: The same as in Fig. but for tan/5 = 4. 



widths, one has a rather large value of the radiative decay width, due to the presence of light 
charginos (low M2) in the W^/x"^ loops, irrespective of the sfermion masses. For instance, 
the point Mi = 180 GeV, M2 = 60 GeV in Fig. 7(a) corresponds to a total X2 decay width 
of about 100 eV and only about 5% of it is due to tree-level decays (in particular, cascade 
decays with m^o = 97 GeV, m-± = 93 GeV, and m^o = 85 GeV). In contrast, in the "unified" 
point Ml = 85.2 GeV, M2 = 170 GeV, one has rTOT(X2) - ^ixl ^ xll) ~ 1 eV for 1 TeV 
sfermion masses, while the total X2 width approaches 50 eV for the lighter sfermion masses 
used in Fig. 7(b), about 90% of it coming from tree-level light-slepton exchange channels. 

A final remark on Fig. 7 is that there are regions where neither the kinematical nor 
the dynamical BR(x2 ^ Xil) enhancement are fully effective, since the regimes for the 
SuSy parameters we outlined in Sees. 3 and 4 are only realized with a large approximation. 
Nevertheless, the combined effect of the two mechanisms can still give rise to large values of 

BR(x^ ^ xh)- 

When tan/? rises (for instance, tan/3 = 4 is assumed in Fig. 8) the physical purity of 
neutralinos decreases (see the corresponding Fig. 4), and mainly effects connected to the 
kinematical enhancement in the scenario (|T^ survive. This can be easily verified by com- 
paring Fig. 8 with Fig. 6(b). The value of the sfermion masses is then a relevant parameter. 
Assuming lighter scalars [Fig. 8(b)] considerably reduces the region of large radiative BR with 
respect to the regime where the sfermion exchange in the tree-level X2 decay is suppressed 
by heavy scalars [Fig. 8(a)]. 

We consider now the effect of varying the parameter fi. From now on, we will keep 
the sfermion masses rather heavy [i.e., we set mg = rrig = 1 TeV] in order to optimize the 
BR(x2 Xil) enhancement. In Fig. 9, the effect of decreasing /i down to ~2Mz is shown 
for tan/3 = 1.2. Again, some insight of the BR behavior can be gained by looking back at 
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Figure 9: Contour plot for the neutralino radiative decay BR in the case tan (3 = 1.2 and 
fi = —2Mz ■ The sfermion masses are all taken degenerate and equal to 1 TeV and = 300 
GeV. The shaded region is excluded by LEPl-1.5 data. 



Fig. 2 (physical "purity" of the neutralinos) and Fig. 5(a) (neutralino mass degeneracy). The 
area around Mi ~ M2 — —fi is again particularly promising. However, quite large radiative 
BR's can be also achieved in three strips of the (Mi,M2) plane, corresponding to different 
kinematical enhancement scenarios [cf. Fig. 5(a)]. Some effects from the reduced dynamical 
mechanism with heavy sfermion masses are also present. 

Keeping fi in the negative range, we now go to larger fi values, which are of particular 
interest for explaining the e'^e~'j'j+ event at the TeVatron. In Figs. 10(a) and 10(b), 
the results for fi = —55 GeV at tan/3 = 1.2 and 2, respectively, are shown. While a large 
BR(x2 Xil) can still be obtained comfortably, we can observe that for increasing tan/? 
the regions excluded by present data extend further and further. 

Going to fi = —Mz/2 at tan/5 = 1.2 (cf. Fig. 11) has the effect of a moderate shifting of 
the large BR regions down to smaller Mi, M2 with respect to Fig. 10(a), accompanied again 
by a drastic reduction of the parameter space allowed by the experimental data. 



6 Top-squark sector influence on BR(x2 XiT) 

As long as m-± ^ Mw, the main contributions to the X2 radiative decay width come from 
the W^/x"^ loops (cf. Fig. 1), unless one considers scenarios with light top-squarks, i.e. 60 
GeV ^ ^ nif (throughout this paper we assume vit = 175 GeV). In the latter case, the 
amplitudes corresponding to t/t loops are non- negligible and can interfere destructively with 
the W^/x"^ loops, hence decreasing the radiative decay width. If one assumes the mass of 
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Figure 10: The same as in Fig. 9, but for fi = —55 GeV and tan/3 = 1.2 (a) and 2 (b). 



the heavy top-squark t2 sufficiently larger than m^, then the bulk of this effect comes from 
the light top-squark ti. Note that the particular values of the other sfermion masses have 
an important influence just on the tree- level neutralino decay widths, hence affecting only 
the BR of the radiative decay. The opposite happens for the top-squark parameters, which 
directly influence the radiative width, and are not involved in the tree-level decays in the 
parameter ranges considered here. In the previous sections, we assumed the two top-squarks 
ti and t2 degenerate, with a mass equal to all the other squark masses. Here, we relax 
this simplification and consider the effects of assuming different values for m^^ and m^^ and, 
in particular, the possibility of a rather light ti. In addition, since the mass eigenstates 
ti 2 can be superpositions of the interaction eigenstates i^^R, we also take into account the 
effect of varying the top-squark mixing angle 9t, defined by: ii = cosOftL + smOttR, t2 = 
— sin Otti + cos 9ttR. 

In Fig. 12, we show the contour plot of BR(x2 Xil) ^ function of the two top- 
squark masses, for four values of the mixing angle 6t and for two interesting choices of 
the neutralino sector parameters which correspond to large-BR(x2 ~^ Xil) scenarios with 
different characteristics. Here, we assume degenerate slepton masses of 300 GeV and the 
other squark masses all equal to 400 GeV. The pseudoscalar Higgs mass is fixed at 300 GeV. 

Here we did not put any restriction on the {m^_^,m^^) plane (besides a rough direct limit 
of about 60 GeV on the top-squark masses from LEP data) and we considered the physical 
masses m^. as independent parameters. However, in the framework of the MSSM, even 
without precise unification assumptions in the scalar sector, one usually derives the physical 
sfermion masses and the L-R mixing angles from the soft SuSy breaking parameters (i.e.: 



"^ijL! T^UR, ''^dn^ •••5 A, ^fe, •••) as well as from n and tan (3 (cf., e.g., Ref.|27[]). Once 



the values of the squark masses of the first two families are fixed to be roughly degenerate. 
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Figure 11: The same as in Fig. 9, but for fi = —Mz/2 and tan/? = 1.2. 



e.g. at 400 GeV as in Fig. 12, and assuming the mass parameters in the sbottom sector to 
have similar values as well (in particular m^^ = mj^), then it seems unnatural to build a 
coherent model with the heavy top-squark t2 considerably lighter than 400 GeV, especially 
when ti is very light. Similarly, the light top-squark is usually lighter than the other squarks. 
Furthermore, one has to take into account the influence of the top-squark sector on the light 
Higgs mass rrih through the radiative corrections (cf., e.g., Ref. |p7||). Especially for low 
tan/3 and for moderate uia, the presence of a very hght ti {mi^ ^ 70-80 GeV) may induce 
nih to fall below the present experimental limits (i.e. 44 GeV or roughly siv?[j3 — a) x 60 
GeV, from Ref.[|TU|), unless t2 is very heavy (m^^ ^ 1 TeV). The value = 300 GeV we 
used in this section generally makes the latter problem negligible, but the discussion above 
still suggests that some regions of the top-squark mass plane we show in Fig. 12 may not 
correspond to a physically acceptable scenario. Nevertheless, having this in mind. Fig. 12 
and the following one provide useful hints to evaluate the pattern of the top-squark sector 
influence on BR(x2 Xil)- 

In Fig. 12(a), we present a typical case of dynamical enhancement (Mi ~ M2 = 60 GeV, 
H = —Mz, and tan/3 = 1.2), while, in Fig. 12(b), a special case of kinematical enhancement 
is considered. The latter (with Mi = 65.2 GeV, M2 = 130 GeV, /i = -70 GeV, and 
tan (3 = 2) does not optimize the kinematical mechanism, but has two interesting features. 
First, the case (b) satisfies the gaugino mass unification relation Mi = | tan^6wM2. Second, 
the value of tan (3 quite different from 1 gives rise to a sizeable mass splitting between the 
neutralinos (~ 17 GeV) and to rather energetic photons, even after the radiative decay of 
a soft X2- Furthermore, in the case (b) a certain amount of dynamical enhancement is also 
present. 

For each fixed value of the radiative BR, we show how the contours move when varying 
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Figure 12: Contour plot for BR(x2 ~^ Xil) the (m^^jm^^) plane. For each fixed value 
of BR(x2 XiT); four curves corresponding to different choices of the mixing angle dt are 
shown with lines of the same style. Going from the lower curve towards the upper one, one 
has: 

(a) Typical dynamical enhancement scenario: 
Ot = -7r/4, (or 7r/2), +7r/4; 

(b) Special kinematical enhancement scenario with gaugino mass unification: 
Ot = 0, -7r/4, +7r/4, +tt/2. 

The other squark masses are taken degenerate and equal to 400 GeV and the slepton masses 
are set at 300 GeV. Also, itlao = 300 GeV. The big black dot corresponds to assuming the 
top-squark masses equal to the other squark masses, as in the previous sections. 



the mixing angle 6t. In particular, for a given BR(x2 ^ Xil), we present four contours 
corresponding to 6t = 0, — vr/4, 7r/4, and 7r/2. This ordering corresponds to going from the 
lower curve to the upper one with the same line style in Fig. 12(b), while, in Fig. 12(a), the 
lower curve is obtained for 6t = — vr/4, the upper one for 6t = +vr/4, and the curve in the 
middle corresponds to the almost degenerate contours for 6t = and 7r/2. As a reference, we 
also show by a big black dot the scenario of complete squark mass degeneracy (in particular 
^ti2 ~ ~ "^^^ GeV), assumed in the previous analysis. Under that hypothesis, one 
would obtain BR(x2 Xil) — 97.5% in the case (a) and BR(x2 Xil) — 65% in the case 
(b). As anticipated, the radiative BR increases with the top-squark masses, showing a larger 
sensitivity to m^^. Indeed, when the mass splitting (m^^ ~ ""^ii) is sizeable, BR(x2 Xil) 
gets independent of the heaviest mass m^^. 

As for the physical composition of the top-squark mass eigenstates entering the radiative 
decay loops. Fig. 12(b) shows that a lightest top-squark corresponding to a pure in {6t = 7r/2) 
is more effective in reducing the radiative BR in the considered case of kinematical enhance- 
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ment. Indeed, for a fixed m^^, taking ti = gives rise to a larger cancellation between 
the t/ti and the W^/x^ loops than in the ti = case. The mixed cases corresponding to 
9t = ±7r/4 fall in the middle of the two pure cases. In the dynamical enhancement scenario 
of Fig. 12(a), the behavior is somehow opposite. It turns out that taking a pure or iji 
state as the lightest top-squark can not be distinguished in BR(x2 Xil) ; while if ii is the 
symmetric (antisymmetric) combination [ii + \/2, the contributions from t/i loops 

are maximized (minimized) and the same happens to the larger destructive interferences 
with the W^/x"^ loops. This is because in the ) the two neutralinos are an almost 

pure photino and an almost pure Hb, which both couple with the same strength to left- 
or right-type standard fermions and sfermions. As a consequence, the contributions from 
pure t/ti and t/tji loops are almost equal, while, if the mass eigenstate considered is an 
antisymmetric combination, its contribution to the matrix element is negligible. When this 
is the case for ti, and ^2 is heavy, the whole top-squark sector does not contribute and the 
process is dominated by W^/x"^ loops. In the case (b), instead, the Xi 2 composition is more 
involved, with sizeable Z components, which distinguish between /l/l and fufn. 

Comparing Fig. 12(a) and 12(b), one notes that varying the top-squark masses can in- 
fiuence the radiative BR more drastically in the case of the kinematical enhancement (b), 
due to the presence of considerably mixed neutralino states, which increases the relative 
importance of the amplitudes from t/t loops with respect to the ones from W^/x"^ loops. 
For instance, assuming m^^ ~ 1 TeV, increasing m^-^ from about 100 GeV to about 1 TeV 
enhances BR(x2 Xil) by a few percent, in the case (a). On the contrary, in the case (b), 
the same growth of the lightest top-squark mass gives rise to a wide increase (of order 100%) 
of the radiative BR. 
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Figure 13: As in Fig. 12, but for a typical case of kinematical enhancement. In (a) the case 
9t = 0, ti = ti is shown; in (b) 6t = n/2, ti = tR. 
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Such effect is even more dramatic in the case of Fig. 13. Here, a typical scenario with full 
kinematical enhancement is shown, with Mi = 178 GeV, M2 = 250 GeV, /i = —2Mz, and 
tan/5 = 1.2. The corresponding, nearly degenerate, neutralino masses are m^o = 180 GeV 
and m^o = 185 GeV. Compared to the cases shown in Fig. 12, in the scenario of Fig. 13 the 
relative weight of the W^/x"^ and t/t loops in the radiative decay matrix element is rather 
different. Indeed, the W^/x"^ and t/t amplitudes depend on (m^i/M^/)^ and {m^/mtY, 
respectively (cf. Refs.p, |l6l)) and here one has m-± as large as 196 GeV. Thus, the t/i and 
W^/x^ amplitudes turn out to be of the same order of magnitude, with different sign. Then, 
one expects that for m^^ ~ 300 GeV (or possibly more, if the heavier top-squark is also rather 
light) the two contributions tend to cancel each other, drastically reducing BR(x2 Xil) 
at the level of a few percent. On the other hand, when the top-squarks, and in particular ti, 
are much lighter (heavier) than 300 GeV, the t/i (W'^/x^) loops dominate, the destructive 
interferences are negligible, and the BR(x2 X1I) can comfortably approach the 100% 
level. This effect is clearly visible in Fig. 13. The influence of a different top-squark mixing 
can be extracted by comparing Fig. 13(a) and Fig. 13(b), where the pure cases ti = ti and 
^1 = in, respectively, are considered. 

7 Conclusions 

In this paper, we showed that SuSy scenarios where the radiative mode for the next-to- 
lightest neutralino decay is dominant do exist and can be naturally realized, especially when 
relaxing the electroweak gaugino mass unification assumptions at the GUT scale. We found 
that very large BR's for X2 ~^ Xi7are obtained when tan /3 ~ 1 and/or Mi ~ M2, with nega- 
tive /i. When Ml = | tan^^vv^2, it is still possible to achieve a large radiative BR, provided 
tan/? ~ 1 and /i < 0. Two different mechanisms, which have different phenomeno logical im- 
plications, can be responsible for the radiative BR enhancement. The dynamical mechanism 
may give rise to signatures including a hard photon and missing energy at hadron colliders, 
where X2 can be copiously produced either in association with a xt/x'j or through the de- 
cays of heavier sfermions. On the other hand, at LEP2, in a scenario with BR(x2 ~* X1I) 
dynamical enhancement, the main neutralino production channel e+e" X1X2 (giving rise 
to a 7-t- ^ signature) is quite depleted, requiring either gaugino or Higgsino sizeable compo- 
nents in both Xi and X2- larger cm. energy e+e" colliders [like the proposed Next Linear 
Collider (NLC)], the dynamical enhancement of the radiative X2 decay may be relevant when 
the X2 is produced in association with a x'j with j > 1, giving rise to signatures containing 
at least a hard photon and missing energy. On the other hand, the kinematical enhancement 
(which implies rather degenerate light neutralinos) can give rise to sizeable rates for final 
photons of moderate energy in both hadron and e~^e~ collisions, provided the next-to-lightest 
neutralino is produced with energy large enough to boost the rather soft photon. 

Finally, an enhanced X2 decay into a photon might be of relevance for explaining, in the 
MSSM framework, the e~^e~'~f'j+ event observed by CDF at the Fermilab TeVatron, that 
is presently one of the most interesting hints for possible physics beyond the SM. 
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